Efficient quantum computation of high harmonics of the Liouville density distribution 
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We show explicitly that high harmonics of the classical LiouviUe density distribution in the chaotic 
regime can be obtained efficiently on a quantum computer ||l]Q. As was stated in Q, this provides 
information unaccessible for classical computer simulations, and replies to the questions raised in 
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In our Letter ^ we showed on the example of Arnold 
cat map that classical chaotic dynamics of exponentially 
many orbits can be simulated in polynomial time on a 
quantum computer. The Liouville density distribution 
P{xi,yj) is encoded on a discretized lattice (2"'' x 2"«) 
using Suq — 1 qubits organized in three registers. After 
each map iteration, the distribution is coded in the quan- 
tum state XlijOylcCi > luj > |0 > with P{xi,yj) = aij 
and Xi = i/N, yj = j/N, N — 2"« ||]. One measurement 
of qubits in this basis gives one point in the phase space 
and therefore the distribution P{xi,yj) can be obtained 
approximately in polynomial number of measurements. 
However the same information can be obtained via classi- 
cal Monte Carlo simulation with a polynomial number of 
orbits, as it was discussed by us in and later repeated 
in 1^ . Based on this observation, the comment makes 
a general claim that no new information can be extracted 
efficiently from quantum computation of such classical 
maps {paragraph 3 in ^) (see also the comment [Q). 
Here we show that this statement is incorrect. Indeed, 
the quantum Fourier transform (QFT) of aij provides 
nondiagonal observables 0|, namely the Fourier com- 
ponents P{kx,ky) = Y,-jex:p{i2T:{kxXi + kyyj))atj/N. 
They obviously contain important information relevant 
for the classical dynamics, and require 0{n^) operations 
including measurement. On the contrary, all known clas- 
sical algorithms will require exponential number of op- 
erations to obtain correct probabilities at high harmon- 
ics k,j;^y ~ N. Such harmonics are important since due 
to chaos a significant part of total probability is trans- 
fered to wavevectors with k ~ exp{ht) where h is the 
Kolmogorov-Sinai entropy, and t is the number of iter- 
ations (see Q). We note that the claim of Q appHes 
equally to the Shor algorithm, where all information is 
also encoded only in squared moduli of amplitudes, but 
where the QFT produces classically unaccessible infor- 
mation. 

In Fig.l we present the probability P{kx, ky) in Fourier 
space for different times t. We note that a 2-dimensional 
(2d) QFT can be efficiently implemented by application 
of usual QFT to each register consecutively. The results 
show that P is composed of well-pronounced peaks, most 
of which move with time to high wavevectors k. They 



remain stable in presence of noise in the quantum gates 
(e.g. top right panel in Fig.l is unchanged if 1% noise 
is added in each gate). The location and amplitude of 
these peaks can be extracted from a polynomial number 
of measurements of qubits after the 2d QFT. 
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FIG. 1. Fourier coefficients \P{kx,ky)'^ of Liouville distri- 
bution for —N/2 < kx^y < N/2, initial state as in Fig.l of 
[1]. Left column: cat map at f = 3, 5, 7 from top to bot- 
tom for riq — 10. Top right: same ai t — 5, riq — 7. Mid- 
dle right: \P{kx,ky)\^ for perturbed cat map (see text) at 
t = 5, riq — 10. Peaks are shown by circles; maximal circle 
size marks peaks with 1 > \P{kx, ky)\'^ > 0.1, circles twice 
smaller those with 0.1 > \P{kx,ky)\'^ > 0.01, etc... Bottom 
right: coarse-grained image of \P{kx, ky)\^ (proportional to 
grayness) for the data of middle right panel, rif — 4. 

For the Arnold cat map the dynamics in {kx,ky) space 
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is especially simple, given by kx = — ky, ky = 
2ky — kx (mod N). However, generally this dynamics 
is very complicated. To exemplify this, we simulated the 
perturbed cat map y = y + x + x^, x = x + y (mod 1) 
(Fig. 1). It can be iterated in 0{n^) operations on a 
quantum computer using modular multiplication. In this 
case, main peaks can be seen directly for short times, 
while for larger times a polynomial number of measure- 
ments of the first n / qubits |^ gives a coarse-grained im- 
age of \P{kx, ky)\'^, including very high harmonics, which 
are unaccessible to classical computation 0. 

As concerns the issue of errors in the computation, 
raised in Q, it should be stressed that the main aim of 
101 was to compare the errors natural for classical and 
quantum computers. In fact, we showed that the natu- 
ral/minimal (last bit) errors for classical computer grow 
exponentially with number of map iterations while the 
natural errors in quantum gates do not destroy the time 
reversibility (Fig. 1 of As a further example. Fig. 
1 of 1^ clearly shows that the errors of relative precision 
10^^ (comparable with ordinary precision on Pentium 
III) completely destroy the reversibility of classical dy- 
namics. At the same time, the quantum errors of relative 
precision 10^^ in operations of quantum gates preserve 
the reversibility. Of course, since the quantum algorithm 
simulates the classical discretized map exactly, the last 
bit errors made on a quantum computer lead to expo- 
nential divergence of nearby trajectories and exponential 
drop of fidelity. This is clearly illustrated by Fig. 3 in 
. Thus exponential instability of classical chaos is pre- 
served in quantum simulations. The quantum computer 
has enormous capacity in memory and precision which 
grows exponentially with the number of qubits; thus in 
the quantum case the last bit errors are much less im- 
portant than for the classical computer with its limited 
memory space. However, the quantum computer has its 
own natural errors related to imperfections in gate op- 
erations. The results presented in [|l],||J^ show that the 
quantum simulation is stable with respect to its natu- 
ral errors in quantum gates operations while the classical 
computation is unstable with respect to its natural errors 
in the last bit of dynamical variables. 
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pie a square of size N/2x N/2 can be built by using 2nq — 2 
Hadamard gates acting on the last Uq — 1 qubits of each 
register x,y. 

A. Liclitenberg and M. Lieberman, Regular and Chaotic 
Dynamics, Springer, N.Y. (1992). 
[7] Contrary to what was argued in pj , high harmonics cannot 
be found from the fine structure of "a few closely spaced 
final points" , since the Fourier coefficients involve the full 
distribution. 
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An initial distribution with exponentially many orbits can 
be prepared efficiently on a quantum computer; for exam- 
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